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The spin-dependent interactions in mesons are considered in detail in the frame- 
' work of the Field Correlator Method. Analytic expressions for the spin-dependent 

o ; 

■ potentials in heavy and light quarkonia are derived with the QCD string moment 

> . 

O . of inertia taken into account. Recent lattice data are analysed using these formu- 



lae and the data are shown to be consistent with very small values of the gluonic 
correlation length, < 0.1 fm. The Field Correlator and the Eichten-Feinberg defini- 



r~| . tions of the spin-dependent potentials in the lowest, Gaussian approximation for the 

QCD vacuum are compared to one another and the two approaches are shown to be 



equivalent in the limit of a vanishing vacuum correlation length, whereas for finite 
values of the latter the difference between these two approaches can be explained by 



I the contribution of the higher-order field corelators, starting from the quartic one. 

in 
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g ■ I. INTRODUCTION 

> 

^ ! Precise knowledge of spin-dependent quark-antiquark interactions can be considered as 

■ one of the most important goals in the QCD spectroscopy. It could help one to identify newly 
discovered orbital excitations, calculate hadronic shifts of mesons with better accuracy, to 
understand effects responsible for a suppression of spin-orbital splittings in higher-mass 
mesons and baryons. Nevertheless, up to now some important features of the spin- dependent 
interactions in mesons remain unclear. In first turn it refers to the order of levels inside nP 
multiplets of heavy-light mesons and to small magnitudes of the fine structure splittings of 
higher excitations. 

It was assumed long ago that in heavy-light mesons the so-called inversion of the P-wave 
levels could take place [ij]. For example, in the D- and 5-mesons, the mass M2 of the 2"*" 
(^X3/2p^^ state might be smaller than the mass of the 0^ state. Here, in the notation 
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n^Pj, the symbol j refers to the total momentum of the light quark, j = whereas n, I, 

and J stand for the radial excitation number, the angular momentum (/ = 1), and the total 
meson spin, respectively. The possibility of such an inversion was checked by calculations in 
different models — see, for example, Refs. . Recently the absence of the spin-orbital 
inversions in heavy-light mesons was studied in Ref. j4| where, in a relativistic potential 
model, the chiral radiative corrections were taken into account and it was shown that such 
corrections provided the absence of the spin-orbital inversion. 

At present the most clean experimental data can be extracted from the D- and 5-mesons 
where the order of levels seems not to be distorted by large hadronic shifts. All four members 
of the D{1P) muhiplet are known D*{2350), Di(2422), Di(2427), and D*{2A60). For the 
Dsj{lP) mesons two of them, D*q(2317) and Dsi(2460), have large hadronic shifts, > 100 
MeV jsl, due to the strong coupling to the DK and D*K thresholds 0, [sl, Si- Recently 
Bi, B2 and Bgi^ B*2 have been reported by the CDF and D0 collaborations [lO|], for which 
Ml < M2 (the notation Mj is used for the mass of the meson with the spin J; spin-orbital 
inversion implies that M2 < Mq). Thus in the observed heavy-light mesons the order of IPj 
levels is similar to that in heavy quarkonia, where Mq < Mi < M2, and coincides with the 
order specific for the pure gluon-exchange interaction. In case of light and i^'-mesons their 
experimental masses also satisfy the condition Mq < M2. 

If one defines the ratio = aso/t, where the matrix elements a so and t are the spin- 
orbital and tensor splittings then, for the Coulomb-type potential, C,c = 1-5. Nonpertur- 
bative Thomas precession leads to a suppression of this ratio so that, in heavy quarkonia, 
^QQ ~ l-Oj a-^d it is even smaller for the D-mesons, ~ 0.8. More specifically, 

a5o(exp) ~ 1.04t(exp) (IP, 66), 

aso{ew) ^ 1.03t(exp) (2P,66), 

a5o(exp) ~ 0.86t(exp) (IP, cc), 

a5o(exp) ~ 0.77t(exp) (IP, D — meson), 

that is, in all cases, the spin-orbital splitting is positive and it is either equal or not essentially 
smaller than the tensor splitting: 

< a5o(exp) < t(exp). (1) 
In other words, the spin-orbital inversion is not observed in experiment. This experimen- 



3 



tal fact contradicts the naive estimates one can make in the framework of potential quark 



models. Indeed, for spin-dependent interactions in heavy guar 



the standard Eichten-Feinberg decomposition is valid 11 



12 



ionia (up to the order 1/m^) 
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=0-4) contains both perturbative (P) and nonperturbative 
The static interquark potential Vo(r), together 



NP/ 



where each potential Ki(r) {n 
(NP) contributions: Vn{r) = V^{r)+V^ 
with the potentials Vi(r) and V2(r), satisfies the Gromes relation 14l |. 

V^{r) + Vlir) - V^{r) = 0. 



(3) 



Notice that this relation refers both to the perturbative and nonperturbative parts of the 
potentials V"„(r) (n = 0, 1, 2). 

In the Eichten-Feinberg representation ([2]) the quark masses rhi and m2 actually mean 
"constituent masses" which are not strictly defined. For example, even for the 6-quark, 
rfib can vary from the value ~ 4.5 GeV up to 5.2 GeV in different potential models, while 
the pole mass of the 6-quark, mfe(2-loop) = 4.79(8) GeV is known with the accuracy ~ 80 
MeV 1^. The situation is similar or even worse for lighter quarks: the pole masses are 
mu,d{2 GeV) ^3^8 MeV, m,(2 GeV) ^ 100(20) MeV, ^ 1.40(7) GeV [51, while the 
"constituent quark masses" in potential models take quite different values, for example, 
^u,d — 0.3 GeV. As will be shown below (see also Ref. |l2|), well-defined average quark 
energies fii appear instead of in the generalisation of Eq. ([2]). 

Then, for a heavy-light meson with mi = ruq and m2 = mq ^ rUg, the splittings aso 
and t are usually presented in the form suggested in Ref. [isj]: 



aso 



4m^ 



A + t, t = — ^T, (4) 
rrigmQ 

where in aso the contribution proportional to rrig^ is neglected. For the "linear -|- Coulomb" 
potential, 

Vo{r) = ar - ^— + const, (5) 

and for the bound state with the angular momentum / and the radial quantum number 
n, the factor A in the spin-orbital splitting (jlj) is therefore given by the matrix element 
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{V({r) = —0" in this case — see Eq. ( I33i) below) 

Ani = {-[Voir) + 2V;{r)])ni = Us{r~')ni - cx{r-')ni, (6) 

which can become negative due to the second (Thomas precession) term in Eq. ([6]). Then, 
since the tensor sphtting t is inversely proportional to the heavy-quark mass rfiQ and is 
therefore expected to be suppressed, one could naively expect that aso < and the inversion 
of levels, M2 < Mq, has to take place. In other words, the spin-orbital inversion could 
appear due to a strong negative spin-orbital interaction which stems from the confining 
linear potential. 

More specifically, one can express the positions of the 0"*" and 2"*" levels via the matrix 



elements aso and t as follows Iq : 



M(0+) = Mo - 2a5o - t, M(2+) = Mq + a^o - O.lt, (7) 

and hence the inversion of levels is possible if 

aso < 0, \aso\ > 0.3t. (8) 

However, as was discussed before, the inequalities ([H]) do not take place for the known heavy 
and heavy-light mesons. 

As will be shown below, there are several reasons for that. First of all, in Eq. ([6]), the 
coefficient A appears to be negative and relatively small even for rather large values of the 
strong coupling constant as- Second, since the first term in aso is proportional to fn^^, then, 
as was already mentioned before, the correct definition of the "mass" rhg for the light and s 
quarks is crucially important. A model-independent definition of rhg as the average kinetic 
energy of the quark is given in Refs. [3] • It will be also given below. Furthermore, the 
very mass factors in the spin-orbital terms in Eq. ([2]) are modified, if rotation of the string 
connecting the quarks and the antiquark is taken into account. This type of corrections 
gives a sizable effect in light and heavy-light mesons (in heavy quarkonia the influence of 
string moment of inertia is much less important) [19] and suppresses aso as a whole. As a 
result, aso appears to be positive and comparable with the tensor splitting. 

In other words, in order to explain experimental data one needs to reconsider spin- 
dependent potentials in detail. This is the aim of the present paper. We employ the 



Field Correlator Method (FCM) (20|. First, we use recent lattice data for spin-dependent 
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potentials in heavy quarkonia 2l| and extract the value of the vacuum correlation length Tg 
2o| . We find that the data are compatible with extremely small values of this correlation 
length, less than 0.1 fm 221. This value appears to be even smaller than that estimated 
previously on the lattice 23|. Notice that the physical role of Tg for the phenomenology 
of hadrons is quite important: in particular, for hadrons of the spatial size R and the 
temporal size Tg, the QCD sum rule method can be applied if R,Tq -C Tg, while potential- 
type approaches are valid in the opposite limit, R,Tq ^ Tg. With such a small value of 
the vacuum correlation length we therefore justify the use of potential-type approaches to 
quarkonia, in particular, the approach of the QCD string with quarks at the ends, which 
will be used then in this paper. Furthermore, the moment of inertia of the string between 
the quark and the antiquark is taken into account, and we show how this string moment 
of inertia affects spin-orbital splittings — it appears as an addendum to the quark energy 
in the denominators in Eq. (jlj) and thus it leads to the above mentioned suppression of the 
term proportional to A. 

For the sake of generality, wherever it is relevant, we keep the contributions of the mag- 
netic and electric correlators separately which allows us to comment on briefly the deconflne- 
ment phase of QCD, where the conflning electric correlators vanish and only the magnetic 
ones survive. 

The paper is organised as follows. In Section [TTl we discuss spin-dependent interactions in 
heavy quarkonia and derive the corresponding potentials in the framework of the FCM. We 
also discuss in detail the difference between the FCM and the Eichten-Feinberg deflnitions 
of the spin-dependent potentials in heavy quarkonia. In particular, we demonstrate that 
this difference is due to the contribution of the higher-order fleld correlators and that it is 
suppressed in the limit of the small vacuum correlation length. In Section UTI] we analyse the 
lattice data for the spin-dependent potentials in heavy quarkonia and for the fleld strength 
Gaussian correlators and extract the vacuum correlation length. We demonstrate the latter 
to be indeed small. Then, in Section IIVI we generalise the form of the spin-dependent 
potentials ([2]) for the case of light quarks deflning the mass parameters rhi through the 
averaged kinetic energies of the quarks and taking into account the QCD string moment of 
inertia. We conclude and discuss the results in Section |Vl 
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II. SPIN-DEPENDENT POTENTIALS IN QUARKONIA 



A. Spin-dependent potentials in the FCM 

We start with the meson wave function in the in- and out-states 

m 

written in a gauge-invariant way with the help of the parallel transporter, 



(9) 



(10) 



The matrix F stands for the vertex function and provides the correct quantum numbers of 
the meson. Then the Green's function can be constructed as 



G,, = {¥;^'\x2,y2\A)^':r{x,,yM)),,A 

= (TrS,(x2, x^\A)T^x^,y^)Sgiy,,y2\A)T^^y2, X2))a, 

where Sg and Sg are the propagators of the quark and the antiquark, respectively, in the 
background gluonic field and we discarded the disconnected contribution. One can proceed 



,(m)t. 



then by using the Fock-Feynman-Schwinger representation 
propagators, thus arriving at 



25 



26 



271] for the single-quark 



dsi / ds2DziDz2 e"^^-^^(Tr(mi - Z)i)F(m2 - Dl)V^W{C)), (11) 
Jo 



where Ki^2 are the quark and antiquark kinetic energy terms. 



dn, K2 



S2 



T^Q + -^zI{t2) 



dro 



(12) 



and the spin-independent interquark interaction is described in terms of the Wilson loop 
W{C), with the contour C running over the quark trajectories: 



W{C) = ex^ig j ds^^{z)F^ 



(13) 



with Sfj^i, being the surface element. In Eq. f[T21) . rUg and ttlq are the conventional pole 
masses of the quark and the antiquark, whereas Si^2 and ri 2 are the proper time variables, 
introduced by Fock and Schwinger |25|; zi^{ti) and -22/^(^2) are the paths of the quark and 
antiquark, respectively. 
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The pre-exponential term in Eq. (fTTi) stems from the fermionic Dirac projectors. It can 
be expressed through the derivatives with respect to the surface element and thus it can be 
pulled out from the average over the background gluonic field. This term contributes to the 
spin-dependent interquark interactions 261] (see Appendices lAl and iBl for the details). The 
most economical way to include the spin-dependent terms in the exponent in Eq. (|TT1) is to 
extend the differential in Eq. (fT3!) to include the spinor structure [27|: 



(14) 



with a 



fjLV 



4i (^M^*^ 



1,2 for the quark and antiquark, respectively). Then the 



spin-dependent interactions result from the mixed terms containing the spin variables. 

The averaged Wilson loop (fT3|) can be expressed through the correlators of the field 
strength tensors as 



exp 



(Triy(C)) 



n=l 



(Tiexpig J dTT^^{z)F^^{z)) 
dTTil)... [ rf7r(n)((F(l)...F(n))), 



(15) 



where the cluster expansion theorem was used (see Ref. 20] for the relevant references 
and for the detailed discussion). The average {{■■■)) stands for connected correlators, for 
example, for the bilocal correlator, ((F(1)F(2))) = (F(1)F(2)) - (F(l)) (F(2)), and F^^ = 
df^A^ — di^A^ — ig[A^, A^] is the vacuum field strength. Obviously, due to the 0(4) rotational 
invariance and colour neutrality of the vacuum, {{F)) = (F) = 0. 

In the Gaussian approximation for the vacuum, when only the lowest, bilocal correlator 
is retained one has, with the accuracy of a few per cent (see Ref. 28|] for the discussion): 



{TrW{C)) oc exp 



dTXf,y{x)d'Kxp{x')Df,yXp{x - x') 



(16) 



where 



Df,uXp{x - x') 



This bilocal correlator of gluonic fie 
scalar functions D{u) and Di{u) as 201 ] 

D^uXpiu) = {S^xSup - Sf,p6i,x)D{u) 



■{{TrF^,{xMx, x')Fxpix'Mx', x))). 
ds can be expressed through only two gauge-invariant 



_d_ 



{ux5up - Up5xu) + 





^11 




pj . 



D,{u). (17) 
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The correlator D{u) = D{uo, \u\) contains only a nonperturbative part and it is responsible 
for the QCD string formation at large interquark separations. The fundamental string 
tension can be expressed as a double integral: 



oo 



a = 2 du dXD{u,X). (18) 
Jo Jo 

The spin-dependent terms in the interquark interaction are generated by the combination 
(^iiuF^y present in Eq. ( |T5l) . which reads: 




(^^..F,, = \ , (19) 



and therefore one needs correlators of the colour-electric and colour-magnetic fields, as well 
as mixed terms, separately. They immediately follow from the general expression (ITTl) and 
read 

|-((TrE,(x)<|.E,(y)$t)) = 5,, [d''{u) + (n) + j + (20) 

|-((Tr/7,(x)$if,(y)$t)) = 5,, [d^{u) + Df(n) + j - (21) 



|^((Trif,(x)<l>E,(y)$t)) = e,,kU,Uk-^. (22) 

where = — y^, u'^ = u^u^. We keep here the superscripts E and H in the correlators 
D and Di in order to distinguish the electric and magnetic parts of the correlators and thus 
to be able to consider a nonzero temperature T. Indeed, while = and Df = = 
Df^ at T = 0, at higher temperatures they behave differently. In particular, above the 
deconfinement temperature, T > Tc, the electric correlator disappears, whereas and 
the magnetic correlators survive. 

The spin-independent interquark interaction is generated by the term oc ds^u{x)dsxp{x') 
in the Wilson loop (|T6l) . It is given by the standard Nambu-Goto action for the minimal 
surface area, 

^string = 0-5'min, 'S'min = / dt j dfi ^/ {ww')^ - wV , (23) 

Jo Jo 

where for the profile function of the string Wfj,{t,P) we adopt the straight-line ansatz: 

w^{t, (3) = (3xi^it) + (1 - /5)x2^, (24) 



9 



Xi 2(t) being the four-coordinates of the quarks at the ends of the string. This approximation 
is vahd at least for not high excitations due to the fact that hybrid excitations responsible 
for the string deformation are decoupled from a meson by the mass gap of order 1 GeV. 

It is straightforward now to derive the spin-dependent interactions in heavy quarkonia in 
the framework of the FCM. For example, the spin-orbital interaction results from the mixed 
terms ds^^axo-dr in (fT6|) and reads: 



so= I dSf,^{w)dTia''^^D^^xp{w - xi) + (1 ^ 2), 



where 



ds 



e daW^it, l3)dbW^{t, /3)dtd/3, a,b = {t, p}. 



For the straight-line string ansatz fl24l) this yields: 

dsii = ridtdp, dsik = eikmPmdtdp, p = [r x {^Xi + (1 - P)x2) 



(25) 



(26) 



(27) 



and thus the angular momentum appears from dsik- Further details of the derivation, as 



well as the comparison with the lattice data, can be found in Refs. |27l . I29l |. The result 
can be presented in the form of Eq. ([2]) with the following identification of the potentials 
{D{u)=D{X,u), Di{u) 



Vlir 
Vi{r 

VAr 



Di(A,z/), u = VA2 + Z/2): 

/»oo nr POO 

2 du dXD^{X,u)+r dvDf{r,v), 
Jo Jo Jo 



-2 



du j d\\l 



- dv j XdXD^iX, iy)+r dvD^ 
Jo Jo Jo 



(r, v) 



(28) 



-2r 



d_ 



duD^{r, u) 



poo 

6 du 
Jo 



1 + -r 



d 



3 9z/2 



Df(r,z/) 



and with the masses fhi replaced by the dynamically generated effective quark masses /ij, 
corrected to account for the inertia of the QCD string — see Section HVl below for a detailed 
discussion. In Appendices |A] and [B] we give the derivation of the terms Vq, V/, and V^' via 
field correlators with the string moment of inertia taken into account. 

With the explicit form of the potentials ( |28l) found in the framework o the FCM one can 



10 



check the Gromes relation ([3]), which now reads: 

du / dX[D^iX,u)-D''iX,u)] 

- ^oo (29) 
+r / di/[Df(r,z/) -Df(r,z/)]. 
Jo 

The right-hand side of Eq. fl30|) vanishes at the temperature T = 0, when D^{u) = D^{u) 
and Df{u) = D^{u) SOj. Notice that, in the FCM, the Gromes relation holds both for the 
perturbative and nonperturbative parts of the potentials ( l28l) . For T 7^ 0, in general case, 
7^ and the Gromes relation does not hold since 0(4) (Lorentz) invariance is violated. 
Therefore lattice measurements of the Gromes relation at finite temperatures, especially in 
the deconfinement phase of QCD, could shed new light on the nature of the vacuum fields. 

The correlator D{u) does not contain perturbative part 31| (note that, for this reason, 
V{{r) — see Eq. f l28l) — does not contain perturbative contributions either) and decreases 
in all directions of the Euclidean space. This decrease is governed by the gluonic vacuum 

has an exponential form, as was found analytically in Ref. 321] and 



correlation length Tg anc 
measured on the lattice 



231]. Then, this correlator has the form: 



D u] 



exp[-M/T, 



(30) 



where the coefficient is chosen to satisfy the relation ( 1181) . The form ( 13011 is oversimplified. 
In particular, it is not regular at m = 0, while a more suitable form should behave as oc at 
small distances \u\ < Tg and then it exponentially decreases at large u's — this follows from 
the gluelump Green's function studied in Ref. [32|. However, this small change with respect 
to the form (l30l) is not felt in the integrals of Vn{r) and cannot be noticed in comparison of 
the field correlators with the lattice simulations, since the region \u\ < Tg is not measured 
on the lattice. The regular Gaussian form was used in Ref. 33|] and both forms were shown 
to yield similar results. 

Another important comment concerning the correlation length Tg is the scale at which 
it is defined (see, for example, Ref. 34| for the discussion of the issue in relation to the 
Operator Product Expansion in QCD). By natural arguments we expect the scale of Tg to 
be of order of the average size of the gluelump (average momentum), which is of the order 
of 1 GeV (see Ref. 

The other correlator, Di{u), contains both perturbative and nonperturbative contribu- 
tions. Its perturbative part leads to the colour Coulomb interaction between quarks, whereas 
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its nonpeturbative part was found through the gluelump Green's function in Ref. [32]: 

D,iu) = D^iu) + Oriu) = ^1 + -^exphVT;], (31) 
where, similarly to Eq. (ITS!) . we have defined 

j'co poo 2 

ai=2 du dXD^^{u,X) = 6TTasa, as = —. (32) 
Jo Jo 47r 

Notice that, for the sake of generality, we consider two correlation lengths, Tg and T^, for the 

correlators D(u) and Df^{u), respectively (see Ref. [3^ for the details) and we find them 

to be different [22 1. 



For large interquark separations, r ^ Tg,Tg, one can find the following asymptotic for- 
mulae for the potentials 



Kir) = ^ + V;{r) = -a, V^{r) = Vs{r) = ^, V,{r) = ^-^nas6^'\r). (33) 

In particular, the static interquark potential comes out from Eq. (I33p in the standard "lin- 
ear + Coulomb" form, 

4 as 

Vo{r) = VQQ{r) =ar- h const. 

O T 

B. Eichten Feinberg versus FCM definition of the potentials V{{r) and V2{r) 

In this chapter we compare the definition of the potentials V{{r) and ^'(r) in the FCM, 
as given in Eq. (128|) . with the Eichten-Feinberg (EF) definition. 



r dr 
VkdV^^ 
r dr 



Eijk^im [ dt t{{g^Hi{n,ti)Ej{r^,t2))), (34) 
lim [ dt t{{g^Hi{ri,ti)Ej{r2,t2))), (35) 



- I) Vrir) + ^fV^ir) = \im ^ dt {{g'H,{r,M)H,{r2M))) . (36) 

with t = t2 — ti. The average is introduced as 

((0(1)0(2))) = Tr(0(l)0(2)W^(5o)), (37) 

with 0(i) (i = 1, 2) meaning the insertion of a field operator (plaquette) in the unperturbed 
rectangular Wilson loop W{So). For the sake of simplicity we consider the limit of the 
antiquark being infinitely heavy, mg = 1712 = oo. 
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Let us start from the quark-antiquark Green's function (ITT!) and notice that the magnetic 



part of the spin-orbital interaction appears from the average (see Refs. 29|, |33|] for the 
details) : 

{gH,{r,,h)W{C)). (38) 

In the last equation, the Wilson loop W{C) is taken along the contour formed by the 
trajectories of the quark and the antiquark, and the averaging is assumed over all such 
trajectories as well as over all vacuum field configurations. For heavy quarks, the deflection 
of the perturbed path C from the straight-line path Cq is small (see Fig. [1]), so that one 
can envisage a sort of perturbative expansion in powers of this small deflection. Below we 
shall perform such an expansion in two ways: 1) in the manner of the Eichten-Feinberg 
approach; 2) in the framework of the FCM. In the given two approaches one uses different 
choices of the surface bounded by the contour C. Indeed, in the FCM this is the minimal 
surface, Spc = 5'min(C'), approximated by the straight-hne ansatz — see Eq. f l24l) — while, 
in the EF approach, the surface is given by the sum of three pieces (see Fig. [1]), 

SEFiC) = So + 6Si + 6S2. (39) 

Due to the nonabelian Stokes theorem, the Wilson loop average ( 138|) does not depend on 
the choice of the surface, hence 

{gH,{n,h)W{SEF)) = {gH,{r,,h)W{SFc))- (40) 
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One can conclude therefore that the two approaches are equivalent, provided the full infinite 
set of field correlators is taken into account. In the meantime, in both approaches, only the 
leading correlators are retained, and this leads to a discrepancy between the definitions of 
the potentials discussed before. Below we explain this in more detail. 

In the FCM one simply uses the cluster expansion for the right-hand side of Eq. (I40l) and, 



keeping only bilocal correlators, arrives at (see Refs. [33[ for the details): 

{gH,{n,t,)W{Spc)) ~ ^ j ds,,{w){g'H,{n,t,)F^,{w)) . (41) 

Then, following the same steps as in Appendix [B], it is easy to arrive at the formula (notice 
that, for the sake of simplicity, we omit here the string inertia and consider heavy quarks): 

which coincides with that given in Eq. (|28|) . This result is accurate up to the quartic 
correlators {{F'^)) and higher. 

Let us now turn to the EF case — the left-hand side of Eq. ( |40|) . Notice that the pieces 
5*0 and 6S2 (see Eq. ( I39ll and Fig. [I]) do not contribute to the average. Indeed, in the sum 

{gH,{r^,t^)W{SEF)) = ig^ Yl / ds^,{H,{r,,t^)F^,{w)W{SEF)) 

= i ds^^{g^Hi{ri,ti)Ff,^{w))+t ds^^{g^Hi{ri,ti)F^^{w)) (43) 

Sn J SS9. 



+1 [ ds^,{g^Hi{n,ti)F^,{w)W{SEF)) 
Jssi 

only the last term contributes, since the integral over Sq is odd with respect to the inversion 
of one of the spatial coordinates, whereas the integral over 6S2 vanishes for ti 3> T^. The 
remaining integral over the surface 6Si can be written as 

{gH,{n,t^)W{SEF)) = i{g''H,{r^,t{) [ ds,,E,{u,t2)W{So)) , (44) 

Jssi 

with 

W{So)=eyipig dsni{w)En{w), (45) 

J So 

and it is easy to check that its cluster expansion starts with the triple correlator with the 
leading correction coming from the fifth-order correlators {{F^)). 
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In order to proceed we write 

z{g^Hi{r,,t^)E^{r^,t2)W{So)) = e,,nrn x ^eijJ^z{g^Hi{r^,h)Ej{ri,t2)W{So)), 

dSji = {Xij{t2) - rij)dt2 ~ Xij{t2){t2 - ti)dt2, 
£ijnrnXlj{t2) = [r X Xii] = — , 

irrii 

where Li is the quark angular momentum. Then 

ALfo = -^J dnaiTT{gH{ruh)W{SEF)) 

/dt 
^^dsj^aiiTT{g^Hi{ri, h)E,{ri, t2)W{So)) 

f / dt, [ dt2it2-h)e,,k'-^TT{g^H,{r,,ti)Ej{n,t2)W{So)) = - [ AVsoir)dti. 

1 J~oo Jtl J 



(TiL 



2m2 

Comparing the last formula with the definition 

^Vsoir) = ^Inir), (46) 

it is straightforward to extract the potential V({r). In particular, after expanding the expo- 
nent in 1^(5*0) and keeping only the lowest contribution, one arrives at 

'fll-}-l = -j tdt J dw^J d\{E{2)E{w)H{l)) + 0{{{E^))), (47) 

where the surface element dsn4 was written as 

T 

dsni = dwndw^i = rndf3dwi = —dXdw^ (48) 

r 

and we defined the triple field correlator, 

{E{2)Eiw)H{l)) ^ e,k,r'^g'{Et{r,M)El{r, - X,w,)Ht{r,M)) ■ (49) 

Notice that each field operator is in fact sandwiched between the two parallel transporters 
(see Eq. like 

E{x) = E{x,xo) = ^{xo,x)E{x)^{x,xo), (50) 

which provide the overall gauge invariance of the field correlator. It is convenient to choose 
the reference point xq = (ti,ri). 

We now employ the relation between the triple field correlator and the derivative of the 



bilocal field corelator [35j: 

-^D{x-y) = ]- r dZna{z){E{x)E{y)H{z)), a{z) = '-^. (51) 
dxn 2 Jy x-y 
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By choosing 

n = 4, t = z - y, y = t, z = t2, x = w^, 
one can rewrite Eq. ( l47l) as 

dr 



^ ^ ' = - tdt / dwi / dX{E{2)E{w)H{l)) 
Jo Jti Jo 

roo Q 

= -2 / (w4-ti)- — D"(w4-ti,X)dwd>^ (52) 



roo 



noo roo pr 

2 D{wi-ti,X)dwid\ = -2 dv I d\D^{u,\). 
Jti Jo Jo 



This result coincides with the FCM expression given by Eq. fH21) . if one neglects the term 
containing A/r. As seen from Eq. fl52l) . this term does not appear in the lowest order of the 
Eichten-Feinberg representation. It is easy to trace the root of this term. In view of the 
equality (HOj) . one can write: 



dVf^ir) ^,,,^5,,, _dVf^{ 



r] 



+0(((F^))) = ^^^^ + 0(((F^))) (53) 

or, keeping only the lowest correction, 

^[Vr{r)-Vr{r)] = 0{{{F''))). (54) 

In other words, the appearance of the correction A/r signifies contributions of higher-order 
terms in the cluster expansion. This correction is clearly of order of Tg/r (see also Eq. (!60|) 
below) and thus it decreases with the decrease of the correlation length Tg. The two rep- 
resentations, the Field Correlator representation and the Eichten-Feinberg one, coincide 
therefore in the string limit of QCD when Tg ^ 0. 

A similar consideration is valid for the potential V2{r), so that one arrives at the following 
formulae for the two potentials in the Eichten-Feinberg representation (to be compared with 
the FCM formulae given by Eq. (EHD): 

dVf^'ir) _ 



dr 
dVj^^ir) 
dr 



poo pr 

2 du dXD^{X,u), (55) 
Jo Jo 

duD^{r,u). (56) 
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Potential 




Vi{r) 






VA{r) 


a, GeV2 


0.22 


0.17 


0.03 




0.07 


as 


0.22 




0.22 


0.16 





TABLE I: The set of parameters, in physical units, extracted from the fits (j57p built for the lattice 



data in Ref. 



2l|. 



III. SPIN-DEPENDENT POTENTIALS AND FIELD CORRELATORS ON THE 

LATTICE 

A. Spin-dependent potentials in heavy quarkonia: The gluonic correlation length 



Recently spin-dependent potentials in a heavy quarkonium were measured on the lattice 
in quenched approximation and without cooling 2l| . Several different lattice configurations 
were used in Ref. 2l|], and the accuracy of these calculations is essentially better than in 



previous lattice calculation^ 
with the results from Ref. 



37 



In this paper we compare the predictions of the FCM 



21| obtained for the lattice size 20^40 and the bare gauge coupling 



13 = 6.0 (this corresponds to the lattice spacing a = 0.093 fm). The fits to the data found 
in Ref. [2li read: 



(T+^, aa^ = 0.0468(2), 



Cv2 



0.297(1), 



0"v2 + 



r 



0.0362(4), 

2 , a,2a' = 0.0070(7), c,2 = 0.288(7), 



(57) 



3c 



v3 



c.3 = 0.214(2), 



/ 2' 
-9 ^g- 



, Crv4 



+ 4—, rUga = 1.16(fixed), g' = 0.292(12), a^^a^ = 0.015(3). 



Parameters of these fits, in physical units, are listed in Table HI One can see that slightly 
different values of the string tension and the strong coupling constant are used in different 
fits in order to reproduce the data. 

With the form of the potentials derived in the FCM and given by analytic expressions 
( 125]) . (|5U]) . and (PTi) we are in a position to compare the predictions of the FCM with the 
lattice data from Ref. (d- In the framework of the FCM. the spm-deperrdent potentials 
are defined with the help of four parameters, {0:5, a, T^, T^}, which possess clear physical 
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as 



a, GeV2 



Tg, fm 



r^, fm 



set 1 
set 2 
set 3 
set 4 
set 5 



0.16 
0.16 
0.16 
0.16 
0.32 



0.22 
0.22 
0.22 

0.17 



0.2 
0.1 

0.07 




0.2 
0.1 
0.1 




TABLE II: The sets of the FCM parameters for the spin-dependent potentials taken from Ref. [21| 
Eqs. ([28|) and ([33j) are used for the sets 1-4 and set 5, respectively. 




r[fm] 



FIG. 2: The profile of the Vq^v) for the set 1 (dash-dotted line), set 2 (dashed line), set 3 (fat solid 
line), set 4 (doted line), and set 5 (thin solid line). Lattice data are given by dots. 



meanings. Notice that none of these four parameters can be treated as a fitting parameter, 
because their values are strongly restricted by phenomenology. Indeed, the string tension 
a defines, for example, the slope of the hadronic Regge trajectories and can be extracted 
from the experimental data for the meson and baryon spectra. The behaviour of the strong 
coupling constant is dictated by QCD, while the values of the correlation lengths are closely 
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related to the gluelump spectrum. Indeed, the inverse of the correlation length gives the 
mass of the lowest one-gluon (Mq^"*) or two-gluon (Mq^"*) gluelump 32]: 



M, 



(1) 



— M^') - - 



(58) 



D^^{u) oc exp (-M(5^Vl) , D{u) oc exp (-M^^^ \u\ 

In the framework of the FCM these masses were calculated to be 2.5 GeV and 1.5 GeV, 
respectively ^] (for the value of the string tension o" = 0.18 GeV^). The relation between 
the gluonic correlation lengths and the gluelump spectrum was also emphasized in Ref. 40 1. 
In particular, using the lattice data from Ref. 4l|, the mass of the lowest gluelump was 
calculated to be 0.9 GeV that corresponds to the gluonic correlation length 0.2 fm, which is 
consistent with the results from Ref. 23||. As will be demonstrated below, the new lattice 



data [21i 
of Ref. 



nve even smaller correlation length, so it would be interesting to repeat the analysis 



40l | for the new data from Ref. 21 1. 
We, therefore, use the FCM form of the spin-dependent potentials, with the parameters 
varying in the range compatible with phenomenology, as a motivated guess and compare 
our theoretical predictions with the lattice data. We consider several such sets of the FCM 
parameters (as listed in Table [Tll), and our aim is to find the set which provides the best 
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0.00 I 1 1 1 1 1 1 1 1 L. 

0.2 0.3 0.4 0.5 0.6 



r[fm] 

FIG. 4: The same as in Fig. [2] but for V2{r). 



overall description of the data. For the sets 1-3 we use the full form of the potentials (1281) : the 
set 4 demonstrates the relevance of nonperturbative interactions since only perturbative part 
of the potentials Vn{r) {n = 0-4) is retained in this case. Finally, the data were approximated 
with the help of the asymptotic large-distance potentials fl5^ (set 5). Comparison of the 
theoretical potentials with the lattice data is given by Figs. In addition, below we also 
compare the lattice data for the field correlators with the predictions of the FCM. 

One can deduce several conclusions from Figs. [2ll6l First of all, by comparing the "purely 
perturbative" set 4 with the rest of the sets, one can see that the data clearly indicates 
the presence of nonperturbative contributions to the potentials, even at small interquark 
separations, where the perturbative physics dominates. Second, by comparing sets 1-3 with 
one another, one can conclude that the data prefer small correlation lengths Tg and Tg. a 
good description is achieved with them both being < 0.1 fm. Then, from the best set 3, the 
gluelump masses can be extracted as Mq^"* = 3 GeV and Mq^^ = 2 GeV, which should be 
confronted with the values 2.8 GeV and 1.7 GeV, respectively, obtained from the predictions 
of Ref. 



32l | after the proper rescaling from a = 0.18 GeV^ used in Ref. [32| to a = 0.22 



GeV^ taken from our set 3. For such small correlation lengths the "large-distance" regime 
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FIG. 5: The same as in Fig. [2] but for Vs{r). The curve for the set 2 coincides with that for the 
set 1. 



of r ^ Tg,Tg happens already for quite small interquark separations, which means that the 
asymptotic large-distance form of the potentials (|28|) should give a good description of the 
data. Indeed, our set 5, based on Eq. (133!) . approximates the data rather well (with the 
exception of the potential V4(r)). For obvious reasons, it fails at small distances. In general, 
the sets 1-3 still give a better description of the data, so one can conclude that the present 
data allow one to study the "anatomy" of the field correlators. 



2l| 



Let us now comment on the Gromes relation ([3]). It was checked numerically in Ref. 
and found to be slightly violated. Indeed, using the data presented in Eq. ( !57|) and in Tabled 
one can find: 

0.018 



Kir) + V;{r) - V^{r) 



\a — a 



vl 



Crv2 + 



Cv2 



0.02 



1 + 



V [fm])" 



GeV^ ^ 0. (59) 



The source of this violation is not clear at the moment. It may stem from an inconsistent 
extraction on the lattice of the potentials entering the Gromes relation which, by itself, 
might be a purely lattice artifact to disappear in the continuum hmit [39|]. Indeed, the 
Gromes relation is exact in the framework of the FCM — see the theoretical expressions for 
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0.6 r 

■ \ 
0.5 - \ 



\ 




-0.1 - r[fm] 

-0.2 L 



FIG. 6: The same as in Fig. [2] but for V4(r). The sets 4 and 5 correspond to delta- functions 
locahsed at r = 0. 



Vn{r) [n = 0, 1,2) given by our Eqs. fl28l) . The Eichten-Feinberg formulae also satisfy this 
relation — see Eqs. (!55l) and (!56l) above for the form of the potentials Vi{r) and V^ir); V^'(r) 
coincides with the FCM one. 

Then, if one adopts the FCM definition (1281) of the potentials V({r) and V^'(r), it is easily 
seen from Fig. [3], that this violation is mostly due to the constant behaviour of the lattice 
potential V({r), which is impossible to get in the FCM formalism with Tg > 0.05 fm. Indeed, 
from the FCM form of the potential Vi(r) given in Eq. (128!) . one finds at r ^ Tg-. 

POO pr / \ \ 2,T'^ AnT 

V;(r) = -2 du d\il--)D{X,u)^-a+^ ^^^'^(0 = -a + (60) 
Jo Jo \ ''^J Jo 

where Eq. fl30l) was used for the profile of the correlator D and only the leading correction in 
the expansion with the small parameter Tg/r was retained, which stems from the term A/r 
under the integral. This correction is responsible for the weak r-dependence of the FCM 
potential Vl{r) depicted in Fig. [31 Notice that a nonconstant behaviour of the potential 
Vl{r) was measured on the lattice in Ref. ^], though a very large uncertainty in the data 
does not allow one to make definite conclusions on the actual r-dependence of this pot ential. 
Consider now the Eichten-Feinberg definitions of the spin-depende nt p otentials [ll| given 



by Eqs. (l34l) - (!36|l and which are used in lattice calculations 2l|, l38l |. The formula ([3 
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implies that corrections, similar to that given in Eq. (l60i) . are not included into the Eichten- 
Femberg defin.fo,. of the potential |V-(r)|'. wh.oh is therefore a constant^ Exactly this 

constant behaviour is confirmed by the lattice calculations [21'] (see the fit V(f^^{r) in the 
set f l57I) ) and, according to Eq. ( l55i) . the potential V({r) is purely nonperturative. On the 
contrary, from Eq. (l56il . one can see that, in the Eichten-Feinberg representation, V2{r) is of a 
purely perturbative nature. In the meantime, a small, though nonvanishing nonperturbative 
contribution to the potential V2(r) was found on the lattice — a nonzero value of the crv2 
in the V^fit('")- Thus we conclude that the problem with the Gromes relation on the lattice 
deserves further investigation with a special attention payed to the surface used in the 
Wilson loop and to the contribution of higher correlators, since the root of the discrepancy 
may reside exactly here, as was explained before. 

Finally, there is a certain contradiction between the theoretical predictions and the lattice 
data for the potential V4(r). Our form of the potential V^lr) given by Eq. (128|) is consistent 
with the delta-functional form of this potential in the limit of the vanishing correlation 
length (see Eq. (1331) ). This should not come as a surprise since the perturbative part of the 
correlator Di{u) reproduces the OGE. Although, for finite values of the correlation length, 
V^lr) is smeared and can become negative, the small-r behaviour of our theoretical curves 
(r < 0.3 fm) is different from that of the lattice data, so this question also deserves additional 
investigation. 



B. Gaussian field correlators on the lattice 



In this chapter we compare the form of the Gaussian field correlators given by the FCM 
with the lattice data taken from Ref. |2l|. As before, we consider the data for the lattice 
20^40 with (3 = 6.0. In particular, we evaluate the correlators of the two plaquettes inserted 
at the positions (0,0,0,0) = (0,0) and (r,0,0,t) = (r,t). Then, with the help of Eqs. (EI]) 
and fl22l) . it is easy to find the correlators measured on the lattice in the form (u^ = + t^): 

^{{9'By{0,0)E,{r,t))) = -rt^^, (61) 
^{{g^B,{0,0)B,{r,t))) = D{u) + Dr{u), (62) 
-^{{9'By{0,0)By{r,t))) = D{u) + D^{u) + r'^^, (63) 
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FIG. 7: The field correlator a'^{{g'^ By{0,0)Ez{r,t))) at r/a = 5 (first plot) and the integral J(r) = 
a'^ {{g'^By{0,0)Ez{r,t)))tdt (second plot). Lattice data are given by dots. 




FIG. 8: The field correlator a^{{g'^B.j:{0,0)B^{r,t))) at r/a = 5 (first plot) and the integral J(r) 
{{g'^Bx{0,0)Bx{r,t)))dt (second plot). Lattice data are given by dots. 



where, in order to comply with the definitions and notations used in Ref. [2l|], we reverted 
the sign of the correlator {{g^ByEz)) (to be compared with Eq. ( l22l) ) and denoted the 
magnetic field as B. Notice that the lattice field strength correlator is calculated as a 
matrix element of an operator defined in (potential) nonrelativistic QCD and simulated on 
the lattice, which is to be multiplied by the appropriate renormalization factor in order 
to give the corresponding correlator in continuum. Although there is an estimate for the 
magnetic factor Zb [42], the electric factor Ze is unknown yet. Since such renormalisation 



24 




factors do not change the shape of the correlators, we adopt the following strategy. The 
profile functions D{u) and Di{u) are taken in the form of Eqs. fl30|) and fl3Tl) . with the set of 
parameters providing the best overall description of the spin-dependent potentials (the set 3 
— see Table ITT] and Figs. W^^- Then the right-hand sides of Eqs. fl^ - flB5|) are maltiplied 
by the normalisation factors A/i, A/2, and A/3, respectively. Notice however that these three 
factors are not independent. Indeed, clearly, J\f2 = A/s since only magnetic fields are involved 
in both correlators. On the other hand, from Eqs. fl6T]) - fl63|) one can derive easily the relation 
between the correlators which is independent of any particular choice of the profile functions 
and fitting constants, namely: 

t[((S,(0,0)5,(r,t)))-((i?,(0,0)i?,(r,t)))] =r((5,(0,0)E,(r,t))). (64) 

This equation establishes the relation between the factors Mi and A/2. In other words, only 
one normalisation factor of the three is independent, for example, A/i. We use it as a free 
parameter. 

The correlators fl6T]) - fl63l) (multiplied by the fitted values of the normalisation factors 
A/i ~ 1.1, A/2 = A/3 ~ 0.6) together with the lattice data are plotted in Figs. [THll From 
these figures one can see a reasonable agreement between the predictions of the FCM and 
the lattice data. The field renormalisation factors can be calculated then as 
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Two comments are in order here. The first comment concerns the small-r behaviour 
of the integrated correlators J(r) (see Figs. 019]), the agreement with the data worsening 
as r —>■ 0. This may be related to two facts. On one hand, the contribution of higher 
correlators at small r's is suppressed by extra powers of r, which might translate into extra 
powers of the ratio r/a on the lattice. Therefore, one should suspect discretisation errors 
of lattice to become significant at short distances. On the other hand, at small interquark 
separations, the perturbative physics dominates, with the behaviour of the correlators being 
essentially dictated by the running strong coupling constant as{r). In our approach adopted 
in this paper we neglect running of the ag which appears to be a good approximation at 
large and moderate distances and fails only at r < 0.2 fm. An accurate description of 
the data at such small distances would require not only considering a running coupling 
constant but also an appropriate modification of the profile functions of the correlators in 
order to warrant a regular behaviour of the latter at the origin, as was discussed before. For 
most applications, however, the large-distance behaviour of the correlators is essential. In 
particular, in the present paper, we apply the FCM to the spin-dependent potentials, where 
also large distances are important. 

The second comment is about the renormalisation factors Ze and Zb estimated above 



and computed on the lattice in Ref. 2l|. Although our colour-magnetic factor Zb ~ 1.3 
appears in a satisfactory agreement with the lattice result Z^* ^ 1.6^1.7, the corresponding 
colour-electric factor Z^^ ^ 1.5 ^ 1.6 differs substantially from the result given in Eq. (|65ll . 
Notice that one should not require exact equality of unobservable quantities, such as the 
field renormalisation factors, measured on the lattice and calculated in the continuum limit, 
because of the discretisation errors and finite volume effects on the lattice. Besides that 
quite different techniques are used in these two approaches, so that establishing a one-to-one 
correspondence between auxiliary quantities may not be possible. Meanwhile, regardless of 
the explicit values taken by the aforementioned renormalisation factors, their ratio is fixed by 
the relation (1641) . This relation is model-independent since it follows from the most general 



parametrisation of the Gaussian field correlator ( |T7I) compatible with its Lorentz structure. 



Then, for the "bare" field correlators calculated in Ref. 



21| to satisfy this relation, the ratio 



las to be Zb/Ze ~ 2 (see Eq. (!65|) ). while it was measured on the lattice to be close to unity 



2l| . This remains an open question which requires further investigation. Concluding this 



discussion let us mention that, since the renormalisation factors for the fields are expected to 
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be calculated directly from QCD, then the relation and similar relations between field 
correlators can be used as a cross-check of the calculated factors or, until these calculations 
are performed, to relate such renormalisation factors with one another. 



IV. QUARK MASSES AND THE EFFECTS OF THE STRING MOMENT OF 

INERTIA 

It was mentioned before that the correct definition of the quark masses is of great impor- 
tance for the correct treatment of the spin-dependent potentials in quarkonia. Our purpose 
here is to define the formal quantities denoted before as mj {i = 1,2). This is especially im- 
portant for light quarks since the very expansion in the inverse powers of the pole quark mass 
is meaningless in this case. In addition we shall take into account the string contribution to 
the canonical angular momentum. 

We return now to the quark-ant iquark Greens's function (JTTil and consider the spin- 
independent part of the meson action. The kinetic terms for the quarks are given by Eq. (fT2l) . 
whereas the spin- independent interaction comes from the Wilson loop fl23|) . 

Using the straight-line string ansatz and considering the system in the laboratory refer- 
ence frame, we synchronise the quark times, 

xio = X20 = t, (66) 

and thus arrive at the standard Lagrangian of the QCD string with quarks at the ends : 



L 



nil V 1 ~ ^1 ^ ""^2 \/ ^ ~ X2 — ar 



[ - [n X iPxi + (1 - P)i2W 

Jo 



(67) 



where r = Xi — X2, n = r/r. The square roots in the Lagrangia n (1671) can be conveniently 
treated with the help of the einbein (auxiliary) field formalism |43|. In this formalism, an 
extra degree of freedom — the einbein field — is introduced. 



mv 1 — a;^ — > 

2/i 2 



yU fix 



(68) 



where the original form of the Lagrangian is restored as soon as the extremum in the einbein 
/i is taken. Thus one reduces the relativistic kinematics to the nonrelativistic one. However, 
the treatment of the system with the help of the einbein field formalism remains a full 
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relativistic treatment since the entire set of relativistic corrections is effectively performed by 
taking extremum in the einbein field. The problem of the centre-of-mass motion separation 
becomes a trivial problem with the einbeins introduced (similarly to the einbeins fii [i = 1, 2) 
for the quarks, a continuous einbein //(/?) is introduced in order to simplify the string term 
in Eq. (|671) ). The result reads (for the sake of brevity we omit the part of the Lagrangian 



responsible for the centre-of-mass motion) 
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4^: 



2/ii 2/i: 







L -- 
where 

R = Cixi + (1 - Ci)x2, r = xi- X2, Ci 



d3 + -u(rnY 



/ii + I3i>d(3 



-/i[r X n]\ (69) 



C2 = l-Ci 



and we have defined the reduced masses for the angular and radial motion separately: 

/il/i2 



(70) 



1^1+ fJ'2' 



(71) 



(72) 



The physical meaning of the variables fii {i = 1,2) is the average kinetic energy of the i-th 



particle in the given state, namely, fii = (A/pM-mf) (see the discussion in Refs. 171, Il8|). 
The continuous einbein variable i^(/3) has the meaning of the QCD string energy density 



17|. 



Following the standard procedure, we build now the canonical momentum as 

p = —— = fiinr)n + fiir — n(nr)), 
or 

with its radial component and transverse component being 



(73) 



(np) = fx(nr), [n x p] = fx[n x r]. 



respectively. Thus we arrive at the spin- independent part of the Hamiltonian 17l |: 

2 



i=l 



fH 

2fii 2 . 



1 

+ / d/3 
'o 



V 

2v ^2 



2/i 2fir 



(74) 



(75) 



One can see therefore that the angular-momentum-dependent term in the Hamiltonian 
5]) contains the total moment of inertia fi which includes both the effective masses of the 
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quarks and also the proper string moment of inertia. Notice that the difference between n 



and fi in the last term in E' 
independent Hamiltonian [1 



5]) gives rise to the so-called string correction in the spin- 



I, 26| and, finally, to the correct Regge slope = 27rcrJ (45 1. 
Clearly, the mass fl appears every time the angular momentum operator is involved. 

We now turn back to the spin-dependent potentials, namely, to the spin-orbital interac- 
tion which must be affected by the string inertia, as was explained before. With the help of 
Eqs. fITOl) and ([71]) one can find that the quantity p, which enters dsik and which gives rise 
to the angular- momentum-dependent interquark interaction — see Eq. (p71) . takes the form: 



p=[rx {(3x, + (1 - f3)x2)] = {13- Ci){r x r] = (/3 - Ci)^. 



(76) 



It is straightforward now to derive the spin-dependent quark-antiquark interactions in 
the full form similar to that given in Eq. ([2]). It is important to stress however that, un- 
like Eq. ([2]), this result is not due to the 1/m expansion but is obtained with the only 
approximation made being the Gaussian approximation for the field correlators. 

Finally, for the spin-dependent interactions, the following modification of Eq. ([2]) can be 
obtained (see Appendices lA] and [B] for the details): 



Ami 



Ami 



(T\Ij <T 2L 



Ami 



Ami 



r dr 



2dVi 
r dr 



1 

2r 



du / dX 



D 



du 



dXD 1 



/il/i 

^) [(1-0)^ 

rj I /ii/i 



X 



((Ti+ CT2)LldV2 1 



2mim2 r dr r' 
{3{crin){(T2n) - (Tia-2) 



dv I XdX 



d 







D + Di + X 



dX' 



C2) 


(T2L 




/i2/(i_ 


C2) 


0-2 -L 








L 


(S) 







,(77) 



12mim2 

(Ti(T2 



(3(crin)(cr2n) - (Ti(T2) 



12mim2 



-VAr) 



6 



dv 



dr 



D(r,u) + 



12/ii/i2 



1 + -r 



d 



3 9z/2 



DAr,u) 



12/ii/i2' 



where the masses fhi are replaced by /ij and fi, which makes this result applicable also to 
light quarks. 

Notice that the potentials Vn{r) themselves do not change their expressions as com- 
pared to Eq. (!28|) . only the "mass" factors are modified. In the limit of heavy quarks, the 
Eichten-Feinberg expression ([2]) is readily reproduced by Eq. fl77j) with the potentials given 
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by Eq. (l28i) . In the meantime, Eq. (177|1 allows one to establish the correction to the inter- 
action ([2]) which comes from the string inertia. Indeed, in the limit rhi ^ (z/), the term jl 
in the denominators can be expanded so that (for equal quark masses, mi = 1712 = m) 



VsDir) 



a 



(78) 



where Vg^j^{r) is given by Eq. ([2]) and S is the total spin of the quark-antiquark pair. The 
term AVsd{^) has the meaning of the string correction for the spin-dependent interquark 
potential, in analogy with the string correction to the spin-independent interaction which 
comes from the similar expansion of the last, angular-momentum-dependent term in the 
Hamiltonian (1751) and which is discussed in detail in the literature — see, for example, 
Refs. [n. 
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46| . The correction AVsnif') constitutes a few MeV for the lowest charmonium 



states. 

The effect of the string inertia is more sizable for light quarks. Indeed, in order to quantify 
the effect in this case let us consider two limits: a heavy-light system, with mi = m, 
m2 — > 00, and the limit of equal quark masses, mi = 1712 = m. In the first limit fi2 ~^ 00, 
Ci ^ and thus the denominator reads (instead of simply /x^): 



/ii 



(79) 



whereas, for the light-light system, one has /ii = /i2 = /io, = 1/2, and the naive term /Xq 
in the denominator is substituted by the expression 



v{(3)d(3 



(80) 



This is exactly the effect we are looking for: the string contributes the total inertia of 
the system together with the quarks and this string contribution is always present in the 
denominator, whenever the angular momentum L appears in the numerator. Due to this 
modification the spin-orbital interaction is weaker than in the "standard" case when the 
moment inertia of the rotating string is neglected. For example, in heavy-light mesons one 
has 

/ii ^ (ar), (81) 

so that, approximating z/(/3) in Eq. fl79|) by (ar), one arrives at the effective denominator 
change 22 1 

1 3 1 
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In the meantime, by an explicit calculation, one can find that, for the P-wave 5-mesons, 
the denominators change constitutes already about 50%, 

1 O.52I. (83) 

Thus we conclude that, for light quarks and at low angular momenta, the suppression of 
the spin-orbit interaction due to the proper string inertia can reach 50%. Notice that this 
suppression has a purely dynamical origin. 

V. DISCUSSION 

In this paper we studied in detail the spin- dependent interactions in quarkonia. Recent 
lattice data for the spin-dependent potentials in heavy quarkonia were analysed using the 
theoretical formulae derived in the framework of the FCM. The comparison with the lattice 
results appears to be in good agreement for the potentials Vq^v), V^'(r), and V3(r). In the 
meantime, for the potentials V{{r) and Vi{r) (at short distances) a discrepancy between our 
theoretical formulae and the lattice data is observed. 

In case of the potential V{{r) this discrepancy comes from the contribution of higher- 
order field correlators and it may be related to the violation of the Gromes relation 
observed on the lattice, since this relation is exact in the FCM for both perturbative and 
nonperturbative parts of the potentials. 

As far as the potential V^lr) is concerned, in the FCM the potential V4(r) stems from the 
delta-function after smearing with a finite vacuum correlation length, whereas on the lattice 
this potential was found to have quite a different shape — it has a "wrong" negative sign as 
compared to its "standard" OGE-inspired behaviour. We conclude therefore that this ques- 
tion deserves a systematic study, primarily on the lattice. If the behaviour of this potential 



reported in Ref. j2l|] is confirmed, this will be a certain challenge for phenomenologists. 

From the comparison of our analytic formulae with the lattice data for the spin- dependent 
potentials in heavy quarkonia we extracted the value of the gluonic correlation length. We 
showed that the data were consistent with extremely small values of the vacuum correlation 
length, less than 0.1 fm. This is an important result since, on one hand, this result ensures a 



very sma. 
vacuum 



1 value of the parameter which governs the cluster expansion in the stochastic QCD 



20|. Indeed, this parameter can be calculated then as 1] = crT^ ^ 0.06 ^ 1 which 
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ensures a fast convergence of the series and justifies working in the Gaussian approximation 
for the field correlators, neglecting higher correlators. On the other hand, such a small 
vacuum correlation length justifies the use of the potential-type approaches to quarkonia up 
to quite small interquark separations. In particular, it validates the use of the QCD string 
approach. Thus we stack to this potential-type approach in order to derive the generalisation 
of the Eichten-Feinberg formula ([2]) for the case of light quarks. In particular, the pole 
quark masses in the denominators are substituted by the averaged quark kinetic energies 
and the string inertia contributes the total inertia of the system and suppresses the angular- 
momentum-dependent interactions. For example, for the spin-orbital interaction, the effect 
of the string inertia may lead to the suppression by almost a factor of two, for low-lying 
states, while, for highly excited states, the spin-orbital interaction is suppressed even more, 
by the factor L"^/^ [2^]. 

In conclusion we return to the quantity discussed in the very beginning of the paper, 
namely to the ratio of the spin-orbit and tensor splittings ^ = aso/t- It was claimed in 
the Introduction that the experimental data on the fine-structure splittings in heavy and 
heavy-light quarkonia appears to be similar, in contrast to the predictions of the heavy quark 
models. The positive sign of ^ and its value ^ ^ 1.0 can be explained if sufficiently large (and 
confirmed by phenomenology) values of the strong coupling constant are used and also the 
correct definition of the dynamical quark masses are introduced. Then, equipped with the 
generalised interquark spin-dependent interactions in the form of Eq. (fTfll . as an example 
we consider the behaviour of the ratio ^ in the heavy-light meson with mi = mq = 5 GeV 
and m2 = mq = 1 MeV. The dependence ^(as) for three realistic values of the string tension 
a is presented in Fig. [TOl From this figure one can see that i) the dependence on the string 
tension is extremely week; ii) for the phenomenologically successful values of the strong 
coupling constant as — 0.5 (see, for example, Ref. [47|), one has ^ ~ 1, in agreement with 
the experimental data, as given by Eq. ([T]). 
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APPENDIX A: COLOUR-ELECTRIC CONTRIBUTION TO THE SPIN ORBIT 

POTENTIAL 

In order to simplify the notations we include the colour factors and the gauge coupling g 
to the definition of the the gluonic field A^. 

We start from the single-quark Green's function and notice that the spin- dependent 
contributions to the interquark interaction come from the Dirac projectors m ~ D and from 
the combination a^yF^y in the exponent (see Eq. (fT5l) ). 

( m-Di io-D \ ( crH crE \ 

m - D = \ , a^yF^y = \ \ = hn + (Al) 

\ -iaD m + D^I \(tE (tH I 

where He = crE, = crH, and for the future convenience we denote the spin-independent 
part of the interaction as ho. Now, in order to build the correction to the diagonal part of 
the Hamiltonian ho + hn which comes from the off-diagonal part /ig, we are to combine the 
latter with the off-diagonal part of the Dirac projector m — D: 



Tr(mi — Di) exp 



- j (iri74(/io + hH + 



(A2) 



Expanding the exponent in Eq. ( ]A2I) and taking the Lorentz trace, one finds that the Hamil- 
tonian has the form: 

H, = ho + hH+^^^^^hE, (A3) 

where the quark energy /ii appeared from D4. Notice that the Hamiltonian Hi defines the 
propagation of the particle in the fifth, Schwinger time ri. Then, after the substitution 
dri = (it/(2/ii), one finally arrives at the following colour-electric contribution to the spin- 
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dependent interaction: 



dx 



Ik 



^ikiCTii [ Xik{Ei{z)Eniu))dSn4iu)— =^ — [ [v X Xi{t)]{E^E^)dtdpdu4 



1-Ci 



/il 



[r X r]n{ExEx)dtdu4d(3 =^ cri 



1-Ci 



1 _ c r°° r' 

{(TiL) ^ / dv d\ 



/ii/i 



r X p]{Ej.Ex)dtdf3du4 



dt = / V^%o{r)dt. 



Supplying the last equation with a similar expression for the second quark, one arrives at 

(A4) 



Vio{r) 



;i-ci)fi^ + (i-c.)p^ 



with 

r dr r jQ 
in agreement with Eq. (1771) . 

For heavy quarks, 



IdV^ 
r dr ' 



(A5) 



1,2), fi^^i 



17111712 



Ci 



mi 



(A6) 



nil + ^1 + 

and thus one arrives at the corresponding term in Eq. ([2]). 

Note that V^'(r) appears only in the first order of this expansion and cannot be exponen- 
tiated — in contrast to V({r) and ^'(r) — see Appendix |Bl This implies that the Gromes 
relations can only be applied to the first order expansion. 



APPENDIX B: COLOR-MAGNETIC CONTRIBUTION TO THE SPIN-ORBIT 

POTENTIAL 

We now turn to the colour-magnetic contributions to the spin-orbital interaction. As 
before, we include the colour factors as well as the gauge coupling to the definition of the 
gluonic field A^. 

Obviously, we consider now the diagonal terms in the matrix (fT9|) . Then one finds from 
Eq. ([15]): 

4o = - / ys'hir)dt= I ds,u{w)dna^^l{Fiu{w)F^p{xi{Ti))) + {I ^ 2) 
^ [ dp [ d\u,\^-^ ainLm{Hr^{w)H^{x{)) + (1 ^ 2), 

2yUl Jo J fl 
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where Eqs. (127|) and (176|) were used and dTi = dt/{2^i), = w^{[3,t) — Xi^(t) {i = 1,2). 
Then, using Eq. fl2T|) for the correlator of the magnetic fields (H^Hn), one can find the 
magnetic contribution to the spin-orbital potential in the form: 



du / dp 



,dDi 



H 



+(1 - Ci)^^ 



,dDl 



(TiL 



(B2) 



2/ii/i 



. + (1-2). 



The two terms in Eq. (lB2p can be combined together with the help of the variable change 
(3 ^ 1 — j3 and, introducing the variable X = (3r, one can arrives at the final formula: 



-— / du I d\ 



r 



+(i-Ci)^^ 



H 



,dDl 



<7iL ^ (T 2-£/ 



(B3) 



which is equivalent to the sum of the potentials Vl{r) and V2{r) — see Eq. 



